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Abstract 

We determine the topological susceptibility y at T = 0 and its 
behaviour at finite T across the deconhning transition in pure SU (2) 
gauge theory. We use an improved topological charge density operator. 

X goes to zero above Tc, but more slowly than in SU{3) gauge theory. 
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1 Introduction 


The anomalous breaking of the flavour singlet axial symmetry in QCD is 
driven by instanton effects |jl|]. This breaking brings about a large mass for 
the pseudoscalar singlet P, P] 

o 2,Nf 9 9 

^v' = + ( 1 ) 

J TT 

X is the topological susceptibility of the pure gauge theory 

X = J d‘^x{0\T{Q{x)Q{0))\0)c^uenched, (2) 

with 

Q(x) = ( 3 ) 

The prediction of eq. (P) is y ~ (180 MeV)^. 

The behaviour of y at high temperature T has also physical relevance (see 
for instance p]). Debye screening in the quark-gluon plasma is expected to 
produce a suppression of the topological susceptibility above the deconhning 
temperature Tc 0. 

In [PI we determined, by a numerical simulation on the lattice, y at zero 
and hnite temperature in the pure SU{3) gauge theory. At T = 0 we ob¬ 
tained the value ixY^^ = 175(5) MeV, which is consistent with previous 
determinations P, §] and with the prediction of eq. i). We also showed 
that y keeps approximately constant below and has a sharp drop at the 
transition point T = Tc- At T/Tc ~ 1.4 the susceptibility y reduces to a few 
per cent of its value before the transition. 

In this paper we determine y for SU (2) pure gauge theory both at T = 0 
and at the deconhning transition with the same technique used in ref. [p| 

In Section 2 we review the method. Our results are presented in Section 3. 
In section 4 we give some concluding remarks. 

2 The method 

The topological charge was measured on the lattice with the operators [y, 

= I: . ( 4 ) 

^ " iiupa=±l 
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^^lvpa is the standard Levi-Civita tensor for positive directions while for neg¬ 
ative ones the relation = —e-^upa holds. IfW jg the plaqnette in the 
fi — h> plane constrncted with i-times smeared links U^\x). These smeared 
links are dehned recnrsively starting from the ordinary link Up{x) as 



The parameter c can be tnned in order to optimize the improvement of the 
operator. We choose c = 0.85 and measnre the topological charge for i =0, 
1 , 2 . 

The topological snsceptibility from the Tsmeared operator on the lattice 
is calcnlated by 

xi” = (EQi‘'WQS;>(o)>. (6) 


X^L niixes with the continnnm snsceptibility y and with all renormalization 
gronp invariant operators of dimension < 4, i.e. the trace of the energy- 
momentnm tensor and the nnity operator Q. The mixing to continnnm y 
is described by the sqnare of the (hnite) mnltiplicative renormalization of the 


topological charge operator to the continnnm operator Q |]l2 


= Z«(/3)Qa" + 0(a®). 


( 7 ) 


Therefore the following relation holds 

xf = Z'>\l}fa(l})\ + + 0 (o«), (8) 

where is the mixing to the trace of the energy-momentnm tensor 

and to the nnity operator. As nsnal (3 = 2Nc/g‘^ where W is the nnmber of 
colonrs and g the gange conpling. 

The additive renormalization can be determined by thermalizing 

the short range flnctnations starting from a zero held conhgnration, with- 
ont changing the (zero) topological content of it |]^ |^, [^, |[. We 
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start with the flat conflguration (all links U^{x) = 1) and create a sample 
of configurations by applying some heat-bath steps. At each step we mea¬ 
sure the topological susceptibility. The content of instantons is checked on 
intermediate steps, by cooling a copy of the conflguration. Configurations 
where instantons or antiinstantons have been produced are eliminated from 
the sample. A plateau is reached after ~ 10 — 20 steps that keeps constant 
along ~ 100 updating steps. The value for is the average of this 

plateau on the ensemble. The number of discarded configurations depends 
on the value of jS. At 13 = 2.5 the rate of discarded configurations per heat¬ 
ing step is ~ 1%; ai (3 = 2.6 it drops to ~ 0.5% and is well below 0.1% at 
(3 = 2.7. 

By subtracting we naturally impose the condition that the phys¬ 

ical topological susceptibility be zero on topologically trivial configurations, 
thus matching the continuum renormalization prescription for y. 

The multiplicative renormalization Z^''\(3) can be determined in a similar 
way. We thermalize short range fluctuations starting from a 1 instanton 
conflguration. At each updating step we measure the topological charge. 
Again checks are performed to eliminate configurations where the topological 
content of the starting conflguration has been changed during the updating 
procedure. is then measured: it stabilizes on a plateau where = 
Z^'^\l3)Q, and the value of Z^^\f3) is extracted from the average on these 
plateaux. 

The statistical accuracy in the determination of the physical value of y 
strongly depends on and Z^^\l3), which in their turn depend on the 

lattice regularization used for the topological charge. The improvement 
of the operator results in a better accuracy Bi- 


3 Results 


The zero temperature determination of y has been done on a 16^ lattice at 
three different values of (3 with the Wilson action. Statistical errors have 
been estimated by using a standard blocking procedure. To fix the scale of 
length we refer to ref. 


17]. We put 


Al 


( 9 ) 
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where f{j3) is the usual 2-loop scaling function and A(/3) is a corrective 
factor tabulated in Table 1 of ref. 0 , which tends to one at large (3, where 
asymptotic scaling sets up. 

The results for x/A^ at zero temperature for 0,1,2-smearings are shown 
in Table I and Figure 1. Scaling is observed in each of them. The three 
determinations (0,1,2-smearings) agree within errors. The smearing process 
has indeed improved the result by drastically reducing the error bars. 

To convert to physical units we need a determination of A/,. We get it 
by using the result of ref. ||^ Tc/A l = 21.45(14), and that of ref. [|^ 
Tc/y/a = 0.62(2) (cr is the string tension). The result is Al = 14.15(42) 
MeV, where, as usual, we have assumed y/a = 440 MeV. Combining the 
2-smeared results at the three f3 values we obtain (y)^/^ = (198±2±6) MeV, 
where the hrst error comes from our determination and the second from the 
uncertainty on Al- 

We have made a comparison with existing determinations of y at T = 0. 
We agree with ref. where (y)^/^=200(15) MeV is obtained by use of 
improved cooling. In ref. 230(30) MeV is quoted for the same quantity, 
which is computed by a different (improved) action and with an improved 
geometric algorithm. The result of refs. 00, when converted in MeV 
by use of the same scale as in this paper, is (178 ± 1 ± 5) MeV. However 
there the renormalization constants were computed by perturbation theory 
because the method of ref. did not exist yet. In ref. y ~ 130 MeV 


is obtained which is lower; we are not able to trace back what part of the 
difference is due to the different method of computing (cooling) and 

what part comes from the determination of the scale a. As for ref. |Q the 
comparison will be done systematically in a forthcoming paper ||^ 


At hnite temperature we used a 32^ x 8 lattice: at this size the deconhning 
transition is located at f3c = 2.5115(40) |^. The temperature T as a function 
of [3 is given by 

The results for 1,2-smearings are shown in Table II and hgure 2: the data 
for the 0-smeared operator have very large errors above the deconhning tem¬ 
perature and are not shown in the hgure. 

At T < Tc our data are consistent with the value at zero temperature, 
while a drop is observed above the deconhning transition. This behaviour 
was also observed in the 517(3) gauge theory H; however in the SU{3) case 
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the drop is quite steeper than in the present case. This qualitative difference 
between the two gauge theories can be well appreciated in Figure 3, where 
x/X(T=o) is plotted versus T/Tc for both SU{2) and SU{3). 


4 Concluding remarks 


We have determined the topological susceptibility of SU{2) pure gauge theory 
at zero temperature and its behaviour through Tc. The improvement of the 
topological charge density operator has made this determination possible. 

If, to £x the scale, the SU{2) string tension is assumed to be the physical 
one, X(T=o) results slightly larger than for SU{3). x is approximately constant 
below Tc, and drops to zero above the transition, however more slowly than 
for SU{3). 
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Figure captions 


Figure 1. x ^.t T = 0. The straight line is the result of the linear ht of the 2- 
smeared data. The improvement from to is clearly visible. 
Ai=14.15 MeV was used to £x the scale by eq. (9). 

Figure 2. x/^i versus T/T^. across the deconhning phase transition for 1 and 
2 smearings. The horizontal band is the determination at T = 0 of 
Figure 1. 

Figure 3. The ratio x/X{t=o) as a function of T/Tc for SU{2) and SU{3) for the 
2-smeared data. 


Table captions 

Table I. x/^i from the 0,1 and 2-smeared operators at T = 0. x*'*^ means the 
continuum susceptibility obtained from the i-smeared operator. 

Table II. T/T„ x^'VAi, X^'VAi as a function of [3. Same notation as in Table I. 



Table I 


p 

10 ^ X 

o 

1 

X 

10 ^ X 

2.44 

4.7(2.1) 

3.73(30) 

3.76(26) 

2.5115 

5.6(2.0) 

3.96(27) 

3.85(19) 

2.57 

3.7(2.1) 

3.97(39) 

3.84(27) 


Table II 


P 

T/Te 

o 

1 

X 

o 

1 

X 

2.40 

0.695 

3.88(46) 

3.45(40) 

2.42 

0.743 

3.42(34) 

3.26(27) 

2.44 

0.793 

3.71(29) 

3.66(24) 

2.48 

0.904 

3.45(28) 

3.53(24) 

2.5115 

1.000 

3.60(31) 

3.54(25) 

2.54 

1.095 

2.99(26) 

2.41(18) 

2.57 

1.203 

2.40(16) 

2.28(12) 

2.60 

1.320 

1.92(17) 

1.81(11) 

2.65 

1.538 

1.46(22) 

1.42(17) 

2.70 

1.786 

1.06(17) 

0.83(8) 
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